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Abstract 

We study the long-range behavior of the heavy quark potential in Coulomb gauge using a 
quenched SU(3) lattice gauge simulation with partial-length Polyakov line correlators. We show 
that the Coulomb heavy quark potential associated with the instantaneous part of gluon propa¬ 
gators in Coulomb gauge, presents a linearly rising behavior at large distances, and the resulting 
Coulomb string tension is greater than the Wilson loop string tension, which can be explained by 
Zwanziger’s inequality. The linearly rising behavior of the Coulomb heavy quark potential persists 
even in the deconfinement phase. The heavy quark potential in Lorentz gauge shows completely 
different behavior than that in Coulomb gauge. Our SU( 3) result, i.e., the Coulomb heavy quark 
potential is confining, qualitatively agrees with that of the SU (2) analysis carried out by Greensite, 
Olejnik and Zwanziger. 
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I. INTRODUCTION 


Recently, the color confinement scenario in Coulomb gauge has been revealing its impor¬ 


tance 


centiy, tne color connn 


Zwanziger [1] discussed the significance of a color-Coulomb 
potential in color confinement. He and his collaborators showed that, in Coulomb gauge, 
the time-time component of gluon propagators, g 2 Doo, including the instantaneous color- 
Coulomb potential and the non-instantaneous vacuum polarization, is invariant under renor¬ 
malization 1,3], where g is a coupling constant of gauge field theory. The instantaneous 
color-Coulomb potential plays an essential role in the Coulomb gauge confinement scenario. 


In the SU( 2) numerical simulation carried out by Cucchieri and Zwanziger 0, it was found 
that the instantaneous color-Coulomb potential D 00 (k) is strongly enhanced at k — 0. More¬ 


over, Zwanziger pointed out that there exists the inequality 


Q 


— Vcoul (R), 


(i) 


where V p h ys (R ) means the physical heavy quark potential and V cou i(R ) the Coulomb heavy 
quark potential corresponding to the instantaneous part of D 00 . This inequality indicates 
that if the physical heavy quark potential is confining, then the Coulomb heavy quark 
potential is also confining. Furthermore, in SU( 2) lattice simulations Greensite et al. found 


that the Coulomb heavy quark potential grows linearly at large quark separations 


6 , Q. 


They showed that the instantaneous part of Doq can be nonperturbatively managed with a 
partial-length Polyakov line (PPL) correlator. (See Ref. Q] for an excellent review.) 

The Coulomb gauge fixing does not fix a gauge completely since it leaves the temporal- 
gauge field free; i.e., one may consider that Coulomb gauge has a remnant residual symmetry 


in the temporal direction. Marinari et al. [fij conjectured that the Coulomb remnant sym¬ 
metry breaking may result in a new order parameter for the deconfinement phase transition. 
Recently, SU(2) gauge-Higgs theory was investigated numerically from this point of view 


a 


In this work, we study the heavy quark potential in the color-singlet channel nonperturba¬ 
tively, using quenched SU( 3) lattice QCD simulations with the PPL correlator in Coulomb 
gauge. We first investigate the behavior of the Coulomb heavy quark potential in the con¬ 
finement and deconfinement phases. The most interesting point is whether the Coulomb 
heavy quark potential has a linearly rising feature at large quark separations. We confirm 
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the consistency between the Coulomb heavy quark potential and the usual Wilson loop po¬ 
tential (or the Polyakov line potential). We repeat the same calculation in Lorentz gauge, 
and compare the results obtained in the two cases. We examine the Coulomb remnant 
symmetry breaking at zero and finite temperatures. 

This paper is organized as follows: In section II, we briefly present definitions of the PPL 
correlators and the Coulomb heavy quark potential in the color-singlet channel. We next 
summarize the gauge fixing methods and the order parameter for the Coulomb remnant 
symmetry. In section III, the numerical results are shown. Concluding remarks are given in 
section IV. 


II. PARTIAL LENGTH POLYAKOV LINE 


In this section, we introduce the Coulomb potential in the singlet channel between two 
static heavy quarks and describe how to fix a gauge on the lattice and the order parameter 
for the Coulomb remnant symmetry. 

Partial length Polyakov lines (PPL) are defined as 

T 

L(x,T)=l[U 0 (x,t), T = 1,2, ■ ■ ■ , N t . (2) 


4=1 


Here U 0 (x,t ) = exp (iagA 0 (x,t)) is an SU( 3) link variable in the temporal direction and a, 
g, A 0 (x,t) and N t represent the lattice cutoff, the gauge coupling, the time component of 
the gauge potential and the temporal lattice size. The PPL correlators in the color-S'f/(3) 
singlet channel are given by 


G(i?,T) = i(Tr[L(H,T)L t (0,T)]), 


( 3 ) 


where R — |x|. From these correlators, we can evaluate the color-singlet potentials on the 
lattice, 

T) = log 

For the smallest temporal lattice extension, i.e., T = 0, we define 



G(R,T) 
G(R, T + a) 


( 4 ) 


V (-R, 0) = — log[£?(/?, a)]. (5) 

The potential V(R, 0) in Coulomb gauge corresponds to a color-Coulomb potential, 
V C oui(R ); the time-time component of the gluon propagator, D 00 (x,t ) = (A 0 (x, i)A 0 (0, t)), 
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can be decomposed into the non-instantaneous vacuum polarization, P(x,t), plus the in¬ 
stantaneous color-Coulomb potential, V cou i(R)5(t), whose term dominates in the Coulomb 
confinement scenario j^|: 


Doo(x,t ) = V cou i(R)5(t) + P(x,t). 


( 6 ) 


In the numerical simulations, the instantaneous contribution has been managed through 
V(R, 0) 0,|?|; this also appears as the enhancement of D 00 at vanishing momentum jsj]. In 
the limit T —> oo, V(R,T) further corresponds to a physical potential, V p h ys , which can 
usually be regarded as the Wilson loop potential in the same limit. In addition, these two 
potentials satisfy Zwanziger’s inequality, Eq. m PI; i.e., if color confinement exists, then 
the color-Coulomb potential is also confining. 

We use the Coulomb gauge realized on the lattice as 


Max EE ReTr u}(x, t) 


(7) 


by repeating the gauge rotations 


2=1 


Ui(x,t ) —>U!f(x,t) — ou^x, t)Ui(x, t)u(x + i, t), 


2.4] ui G SU (3) is a gauge rotation matrix and Ui(x, t ) are spatia 

in 


where 

25I ] Thus, each lattice configuration can be gauge fixed iteratively 


( 8 ) 


lattice link variables. 


The temporal gauge fields possess gauge freedom even after the Coulomb gauge fixed. 
One can still perform a time dependent gauge rotation on the Coulomb gauge fixed links: 


Ui(x,t) —► cJ(t)Ui(x,t)uj(t), 
Uo(x,t) —> u)\t)U 0 (x, t)u(t + 1). 


(9) 


It was conjectured by Marinari et al. |9| that this remnant symmetry in Coulomb gauge 
is closely related to color confinement physics. Recently, Greensite et al. [7] proposed the 
following order parameter for the Coulomb remnant symmetry: 

N t 


1 


Qs 


t =1 


-Tr 


U(t)W(t) 


( 10 ) 


u(t) = 


( 11 ) 
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where V s = N%, and N s stands for the spatial lattice size. The quantity Q s is invariant 
under the transformation ©• Therefore Q s = 0 if the symmetry is not broken. This 
quantity contains 

Y2 V TH 7 „(f, t)U„(y, i)~iy e -v(\*-y\) (12) 

s s 

Consequently, if V increases linearly for large \x — y\, i.e., it is a confining potential, then 
Q s —> 0 as V s —> oo, and if V —> const, as \x — y\ —> oo, then Q s > 0. Hence, we may consider 
that Q s is an order parameter of the confinement and deconfinement phase transitions. If 
this is the case, it would be a very desirable order parameter because this works also for the 
full QCD including dynamical quarks. 


III. SIMULATION RESULTS 

In this study, we carried out SU(3) lattice gauge simulations in the quench approximation 
to calculate the PPL correlators in the color-singlet channel. The lattice configurations were 
generated by the heat-bath Monte Carlo technique with a plaquette Wilson gauge action, 
and we adopted the iterative method [ll| for fixing a gauge. 

A. Coulomb heavy quark potential 


TABLE I: The fitting results for the string tensions. K\ in lattice units and a\ in physical units 
were obtained from the singlet V(R,0). K w and a w stand for the Wilson loop string tension at 


(3 = 6 


.o Q 


(not calculated here). We use the relation = \[Ka 1 and the lattice cutoffs 


estimated with the Monte Carlo renormalization analyses 


(3 IQ [MeV] K w JoZ, [MeV] 

5.85 0.2291(22) 706(4) 

5.90 0.1950(10) 716(4) 

5.95 0.1726(6) 736(3) 

6.00 0.1467(4) 740(3) 0.0513(25) [12] 470(46) 
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V(R,0), 18 4 



( 3=600 

| 3 = 5.90 

( 3 = 5.95 

( 3 = 5.85 


FIG. 1: Coulomb heavy quark potentials obtained from the PPL correlator with T = 1. These data 
were obtained in the 18 4 lattice simulation at (3 = 5.85 — 6.00. We used 400 gauge configurations 
measured every 100 sweeps. 


Figure [T| shows the results of the Coulomb heavy quark potential V(R , 0) obtained from 
the PPL correlator with T — 1. These data were obtained from the 18 4 lattice simulations 
at /3 = 5.85 — 6.00. To obtain the string tension, we assume the following fitting function: 


V(R,T) = C + KR + A/R, A = -nf 12 , 


(13) 


where C is a constant and K corresponds to the string tension. We find that the fittings 
are good; x 2 /nd/ ~ 0(1) for the fitting range R = 2 — 6. It is found that the Coulomb 
heavy quark potential V(R, 0) rises linearly as the distance R increases at (3 = 5.85 — 6.00, 
and hence it can be described by the linear rise function with the string tension. The string 


121 for /3 = 6.0 are listed 


m 


tension for /3 = 5.85 — 6.00 and the Wilson loop string tension 
Table □ Ki at (3 = 6.0 is approximately three times larger than K w , and the value of ^ fa\ 
increases with (3. Although similar results for the 3 dependence of the string tension were 
also obtained in the SU(2) lattice simulations 6,1.7], we have no clear explanation. 


6 








V(R,T), 18 3 x32, P=6 



FIG. 2: T dependence of the Coulomb heavy quark potential at (3 = 6.0. The symbols with error 
bars and the solid curve correspond to the numerical data calculated in the present work and using 
the Wilson-loop potential [r|, respectively. As T increases, V(R,T) seems to approach V W (R). 


Because V(R , T ) in the limit T —> oo corresponds to a physical potential, we expect that 
V(R , T ) becomes comparable with the Wilson loop potential, V W (R), when T becomes large 
enough. The T dependence of the Coulomb heavy quark potential at (3 = 6.0 is displayed 
in Fig. [21 We used the 18 3 x 32 lattice and 600 configurations measured every 100 sweeps. 
V(R, T) may approach V W (R) as T increases. Q The T dependence of the PPL correlator is 
led completely here although there is consistency in the 577(2) lattice calculations 
and [7]. 

The finite temperature behavior of the Coulomb heavy quark potential is shown in Fig. 
13 This simulation was carried out on the 24 3 x 6 lattice at (3 — 6.11 (a ~ 0.086 fin), 
corresponding to T/7) ~ 1.50, where T stands for the system temperature and T c the 
critical temperature of the quark-gluon plasma phase transition. The 300 configurations 
measured every 100 steps are used. It is very remarkable that even in the deconfinement 
region, T/T c 1.50, the Coulomb heavy quark potential V(R , 0) is not screened and still a 


not contro_ 
in Refs. E 
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V(R,T) at finite temperature 



FIG. 3: Linear-rise and screened potentials in the deconfinement phase at T/ T c ~ 1.50. Note that 
V(R, 0) still indicates linearly confining behavior, while the potentials V ( R , T) with finite T are 
screened and are consistent with the calculation employing the usual full Polyakov-line correlator. 


linearly increasing function at large distances. Fitting these data with Eq. m, we obtain 
K\(/3 = 6.11, T/T c ~ 1.50) = 0.118(1), and we find that A fo\((3 = 6.11, T/T c ~ 1.50) = 
792(10)MeV, which is larger than the value y/o\(/3 = 6.0, T ~ 0) = 740(4) MeV in Table 
IU thus, it seems that the Coulomb string tensions with T = 1 depend on QCD coupling 
(3 rather than the system temperature T. In the case of the color SU(2) simulation, it is 


reported in Ref. that the Coulomb string tension scales well according to the two-loop f3 
function. On the other hand, as the temporal extension T increases, the potential V(R,T) 
are screened at large distances, R > 1/T = 6 on this lattice, and they show the usual finite- 
temperature screening dynamics that have been investigated nonperturbatively in lattice 


simulations |14L 1 
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B. Remnant symmetry in Coulomb gauge 


We investigate the role of the remnant symmetry in Coulomb gauge. We calculate the 
order parameter Q s given in Eq. (USD- and here we also study the color average order 
parameter 

<3™ = y E < ■ (14) 

The volume dependence of the Q values at zero temperature is shown in Fig. 01 and the 
simulation parameters are listed in Table UT1 It is found that the 0 values vanish as V s —> oo. 
This indicates that the remnant symmetry is unbroken, and the system is in the confinement 
phase. The Q values in the deconfinement phase are displayed in Fig. 0 The simulation 
parameters are also listed in Table HT1 Figure |3] shows that the Q values at finite temperature, 
T/%. ~ 1.26, also go to zero as V s —> oo. It is found that the Q values calculated here 
indicate no qualitative difference between the confinement and deconfinement phases in 
pure SU( 3) gauge theory; it is surprising that all the Q values vanish as V s —> oo even in the 
deconfinement phase, but this is a consequence of the fact that the Coulomb heavy quark 
potential V (i?, 0) is always confining. 

TABLE II: Simulation parameters for the calculation of Q av and Q s . All configuration measure¬ 
ments were performed every 10 steps. 



Nf x 32 


A 3 x 6 

N s 

No. of conf. N s 

No. of conf. 

8 

100 

18 

100 

18 

100 

24 

90 

24 

100 

42 

30 

32 

30 




C. Heavy quark potential in Lorentz gauge 


Although the argument concerning the Coulomb heavy quark potential defined by the 
PPL correlator f| 7] is based on use of Coulomb gauge, Q we carry out the same calculation 
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N s x 32, [3=6.0, zero temperature 



FIG. 4: Volume dependence of Q in the confinement phase at zero temperature. Q s and Q av 
vanish as V —> oo. 

in Lorentz gauge. The results are displayed in Fig. [HI This calculation was performed on 
the 24 3 x 32 lattice at (3 = 6.0, and 100 configurations measured every 100 steps are used. 

The potential V ( R , 0) in Lorentz gauge represents completely different behavior com¬ 
pared with the case of the Coulomb gauge. It is flat even at large distances. Moreover, 
as T increases, the Lorentz heavy quark potential tends to approach the usual Wilson loop 
potential from below, but at large distances, the Lorentz heavy quark potential seems not to 
be confining. This strongly suggests that the color confinement mechanism is very different 
between Coulomb and Lorentz gauges. 

IV. CONCLUDING REMARKS 

We have nonperturbatively studied the long-range behavior of the Coulomb heavy quark 
potential defined by the partial-length Polyakov line correlators in quenched SU(3) lattice 
gauge simulations. We confirmed numerically that the Coulomb heavy quark potential, 
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N' x 6, [3=6.0, finite temperature 



FIG. 5: Volume dependence of Q in the deconfinement phase at finite temperature, T/T c = 1.26. 
Q s and Q av vanish as V —> oo. 


corresponding to the color-Coulomb instantaneous part of the time-time gluon propagator, 
is confining, as suggested in the Coulomb confinement scenario 1], The SU( 3) results 
obtained in this study are qualitatively consistent with those obtained in the SU(2 ) analysis 
carried out by Greensite et al. jsj, I?]. It is significant that we treated the instantaneous part 
nonperturbatively in the numerical simulations. 

The Coulomb heavy quark potential V ( R , 0) in the confinement phase rises linearly at 
large distances. Its string tension is several times larger than that obtained from the usual 
Wilson loop potential. As the temporal extension, T, of the PPL correlator, increases, the 
Coulomb heavy quark potential asymptotically approaches the usual Wilson loop potential. 
In these simulations, complete agreement is not confirmed. Note that consistency was found 
in the SU( 2) lattice gauge simulation jf|. Furthermore, the result that the SU(3) Coulomb 
heavy quark potential indicates a stronger confining property can be understood from the 
relation V phys (R) < Q- 

The Coulomb string tension can also be estimated from lattice calculations of gluon 


11 








V(R,T) in Coulomb and Lorentz gauges 



FIG. 6: Heavy quark potentials in Coulomb and Lorentz gauges at (3 = 6.0. The symbols with 
error bars represent the numerical data in Coulomb and Lorentz gauges, and the thick solid curve 
represents the Wilson loop potential at (3 = 6.0. 


17] that the Coulomb string tension 

Q. 


propagators in Coulomb gauge. It is reported in Ref. 

&coui ~ cr w n S on■ However, in the extensive lattice study of Langfeld and Moyaerts 
the authors concluded that the result o cou i ~ (2 — 3)a wi i son cannot be ruled out. This 
observation is consistent with results derived in SU(2) and 57/(3) lattice calculations with 
the PPL correlator. 


In 517(2) lattice simulations |7j, (3 dependence of the Coulomb string tension scales as 
a two-loop (3 function. The Coulomb string tension measured here depends on the QCD 
coupling (3 = 6 /g 2 . This (3 dependence is not explicable here, and it is therefore desirable 
to verify it with numerical simulations in higher f3 regions. 

The Coulomb heavy quark potential V(R, 0) in the deconfinement phase at T/T c ~ 
1.50, also increases linearly with R at large distances; i.e., it is not screened. However, 
the Coulomb heavy quark potential with the finite temporal length in Coulomb gauge is 
sufficiently screened, and finally it becomes comparable to the screened potentials obtained 
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from the full temporal length Polyakov line correlator. This result may be explained in 


the Coulomb confinement scenario [1|, |3[. V(R, 0) corresponds to the instantaneous part, 
V cou i(R)m, an d therefore V cou i(R) does not depend on the time (temperature), whereas 
the non-instantaneous vacuum polarization term causes the time (temperature) dependent 
contribution. 

We have also investigated the behavior of an order parameter related to the remnant 
symmetry in Coulomb gauge. As the lattice volume approaches infinity, the order parameter 
Q vanishes in both the confinement and deconfinement phases. Thus, this order parameter 
is not a good parameter for understanding confinement physics, at least in pure SU(3 ) gauge 
theory. This feature is a consequence of the fact that the Coulomb heavy quark potential is 
always confining in the confinement and deconfinement phases. 

The Coulomb gauge plays an essential role in the confinement scenario discussed in Ref. 

1]. We carried out the same calculation for the Lorentz gauge and found that the behavior 
of the Lorentz heavy quark potential is completely different from the case of Coulomb gauge. 
In Lorentz gauge, a different confinement scenario should be called for. 

The color-Coulomb instantaneous part of the time-time gluon propagator exhibits color 


confinement. As discussed in Refs. 


[1] and |3], 


this expectation may be satisfied in a 


dynamical-quark lattice simulation. The vacuum polarization causes a quark-pair creation 
at large quark separations, i.e. “string breaking”, whereas the Coulomb linear-rise potential 
also exists in that case. If one nonperturbatively extracts the contribution of the vacuum 
polarization, it may be easier to see the string breaking. 

The color-Coulomb instantaneous part, V cou i(R)5(t ) or V(R, 0) on the fixed-time slice, 
shows the strong linear-rise behavior in SU(2 ) and SU(3) gauge theories at zero and fi¬ 
nite temperatures; this phenomenon is as expected in view of the Coulomb confinement 


scenario 


io Q, 3- 


In terms of the Faddeev-Popov (FP) operator M, V cou i(R) is given by 
(M -1 (—)M~ 1 ), which has time-independent and long-range properties. According to the 
Gribov picture, gauge configurations for which the eigenvalues of the FP operator are small, 
i.e., regions near the Gribov horizon, play an important role in the confinement dynamics 
Q. Greece et a,, investigated eigenvalue distributions of the FP operator using the 
SU (2) lattice gauge simulation J20j. 

Here we concentrate on the Coulomb gauge. It would be interesting to study how the 
picture changes when we employ the other gauge by interpolating two gauges, such as 
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yif-i diAi + A d t A t — 0 211. 

In this study, we use the gauge rotation, e iadiAi , which is used in the stochastic gauge 


fixing|22|]. This fixing term is “attractive” in the Gribov region and “repulsive” if one 
eigenvalue of the Faddeev-Popov operator is negative, and therefore the effect of Gribov 
copy may not be so serious. Nevertheless, it is important to study this point. 

We investigated the behavior of the color-singlet qq potential in the present calculation. 
However, an extensive numerical study of the color-dependent forces between two quarks may 
be necessary to understand multiquark hadrons. In Coulomb gauge, 577(3) color-dependent 
potentials have also been calculated in lattice gauge simulations|23l]. 

It is interesting from a phenomenological point of view that the linearly rising behavior 
of the Coulomb potential exists even in the deconfinement phase at finite temperature. 
Therefore, a higher temperature simulation is indispensable and it may help to elucidate the 
behavior of the complex quark-gluon plasma system. 
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